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Nonlinear  extensional  vibrations  of  quartz  rods 

H.  F.  Tiersten 

Department  of  Mechanical  Engineering.  Aeronautical  Engineering  and  Mechanics,  Rensselaer  Polytechnic 
Institute,  Troy,  New  York  12181 

A.  Ballato 

U.  S.  Army  Electronics  Technology  and  Devices  Laboratory,  USAERADCOM,  Fort  Monmouth.  New  Jersey 
07703 

(Received  8  October  1982;  accepted  for  publication  4  February  1983) 

The  one-dimensional  scalar  differential  equation  describing  the  extensional  motion  of  thin 
piezoelectric  rods  is  obtained  from  the  general  nonlinear  three-dimensional  description.  Only  the 
elastic  nonlinearities  are  considered.  The  relations  between  the  quadratic  and  cubic  coefficients  of 
the  rod  and  the  fundamental  anisotropic  elastic  constants  of  various  orders  are  derived.  The 
quadratic  rod  coefficients  are  calculated  for  various  orientations  of  quartz  rods,  but  not  the  cubic 
rod  coefficients  because  the  fundamental  elastic  constants  of  fourth  order,  which  are  required  for 
the  calculation,  are  not  presently  known.  The  nonlinear  equation  and  boundary  conditions  are 
applied  in  the  analyses  of  both  intermodulation  and  nonlinear  resonance  of  quartz  rods.  In  each 
instance  a  lumped  parameter  representation  of  the  solution,  which  is  valid  in  the  vicinity  of  a 
resonance,  is  obtained  and  the  influence  of  the  external  circuitry  is  included  in  the  treatment. 

PACS  numbers;  43.40.Cw,  62.20.Dc,  43.40.Ga 


INTRODUCTION 

A  one-dimensional  scalar  nonlinear  differential  equa¬ 
tion  describing  the  extensional  motion  of  thin  piezoelectric 
rods  oriented  in  an  arbitrary  direction  with  respect  to  the 
principal  axes  of  an  arbitrarily  anisotropic  crystal  is  obtained 
from  the  general  nonlinear  three-dimensional  equations  of 
electroelasticity. '  Only  the  elastic  nonlinearities  are  includ¬ 
ed  in  the  description.  The  electrical  behavior  is  taken  to  be 
linear  since  quartz  has  small  piezoelectric  coupling.  The 
treatment  provides  the  relation  between  the  quadratic  and 
cubic  nonlinear  extensional  coefficients  of  the  rod  and  the 
fundamental  anisotropic  elastic  constants  of  second,  third, 
and  fourth  order,  along  with  the  well-known  relation 
between  Young’s  modulus  and  the  fundamental  second-or¬ 
der  elastic  constants.  The  quadratic  rod  coefficients  are  cal¬ 
culated  for  various  orientations  of  quartz  rods  with  respect 
to  the  principal  axes  of  the  quartz  crystal.  Such  calculations 
cannot  be  performed  for  the  cubic  rod  coefficients  because 
the  fourth-order  elastic  constants  of  quartz,  on  which  the 
cubic  coefficients  depend,  are  not  presently  known. 

The  steady-state  solution  of  the  nonlinear  system  is  de¬ 
termined  at  an  intermodulation  frequency^  by  straightfor¬ 
ward  iteration  from  the  linear  solution.  The  intermodulation 
current-voltage  relation,  which  is  valid  in  the  vicinity  of  a 
resonance,  is  obtained  and  the  crystal  is  incorporated  in  the 
test  circuit  and  the  relation  between  the  intermodulation  and 
driving  voltages  is  determined.  This  relation  and  the  calcu¬ 
lated  quadratic  rod  coefficients  can  be  used  along  with  mea¬ 
surements  to  determine  the  cubic  rod  coefficients  for  any 
desired  orientation  of  quartz  rod,  including,  of  course,  the 
practical  S*  X  cut. 

In  addition,  the  nonlinear  equation  and  boundary  con¬ 
ditions  are  applied  in  the  analysis  of  nonlinear  resonance  in 
quartz  rods,  llie  steady-state  solution  is  obtained  by  means 
of  an  asymptotic  iterative  procedure^  and  an  expansion  in 


the  linear  eigensolutions.  A  lumped  parameter  representa¬ 
tion  of  the  nonlinear  solution,  which  is  valid  in  the  vicinity  of 
a  resonance,  is  presented.  The  representation  relates  the  am¬ 
plitude  of  the  rod  displacement  nonlinearly  to  the  voltage 
across  the  electrodes.  The  external  circuitry  is  incorporated 
in  the  analysis  and  an  equation  relating  the  amplitude  of  the 
rod  displacement  nonlinearly  to  the  driving  voluge  and  oth¬ 
er  circuit  parameters  is  obtained. 

It  should  be  carefully  noted  that  the  coupling  between 
extension  and  flexure,  which  can  exist  in  an  anisotropic  rod, 
is  not  included  in  this  description.  However,  since  the  cou¬ 
pling  of  a  high  overtone  of  flexure  with  a  low  extensional 
mode  occurs  only  for  certain  well  separated  geometries,*  the 
equation  for  the  anisotropic  rod  presented  here  is  valid  for  all 
other  geometric  ratios.  This  means  that  the  equation  is  valid 
for  almost  all  geometries  and  certainly  for  all  practical  cases. 

I.  NONLINEAR  EQUATIONS 

The  stress  equations  of  motion  and  charge  equation  of 
electrostatics  for  material  points  of  an  electroelastic  solid 
with  small  piezoelectric  coupling  may  be  written  in  the  re¬ 
spective  forms' 

^U,L  (1) 

I>L.L=0.  (2) 

where  we  have  employed  the  conventions  that  capital  and 
lower  case  indices,  respectively,  refer  to  the  reference  Carte¬ 
sian  coordinates  and  present  Cartesian  coordinates  of  mate¬ 
rial  points.  We  also  employ  the  conventions  that  a  comma 
followed  by  a  capital  index  denotes  partial  differentiation 
with  respect  to  the  known  reference  coordinates,  i.e.,  the 
independent  variables  excluding  time,  a  dot  over  a  variable 
denotes  partial  differentiation  with  respect  to  time,  and  re¬ 
peated  tensor  indices  are  to  be  summed.  The  symbols  p°, 
®nd  D^,  respectively,  denote  the  reference  mass 
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density,  the  mechenicel  displaoement,  the  Fiola-ICifchhoff 
stiese  tensor,  and  the  electfic  displacenient  vector.  In  the 
nonlinear  stress  equations  of  motion.  Eq.(  I),  it  is  understood 
that  the  motkm  of  a  material  point  is  described  by  the  (bnc- 
tional  relation 

(3) 

which  m  ooe-to-coe  and  differentiable  as  often  as  required, 
where and  lespectivdy,  demMe  the  present  and  refer¬ 
ence  poeition  of  material  points.  The  symbol  (1) 

a  translation  operator,  which  serves  to  trai^te  a  vector 
from  the  present  to  the  reference  position  and  vice-versa  and 
is  requited  for  notatkmal  consistency  and  clarity  because  of 
the  use  of  capital  and  lower  case  indi^,  respectively,  to  refer 
to  the  reference  and  presoit  positions  of  material  points. 
Clearly,  the  mechanical  displacement  vector  Uj/  andpreser' 
and  rrference  positions  of  material  points  and  are  re¬ 
lated  by 

For  purposes  of  this  treatment,  the  constitutive  equa¬ 
tions  for  ff|y  andZIf.  may  be  written  in  the  form' 

=yjJt^lM’  (5) 

—  ^UtK  +  <Lt  ^ K> 

where 

^LV  =  ^UHK  ^ro  ~  *KLM  *  +  J  ^IMFOAt  ^FO  ^AB 

+  i  ^U/rOABCD  ^AM  ^CDt  (7) 

and 

^AB  —  ^yi.A  yi,B  ~ ^ab)  (®) 

is  the  finite  strain  tensor.  Substituting  from  Eq.  (4)  into  Eq. 
(8)  we  obtain 

^AB  —  +  **But  +  (^1 

and,  of  course,  since  the  electric  field  is  quasistatk,  we  have’ 

»’s=-f'4c-  (»0) 

In  view  of  the  small  piezoelectric  coupling  in  quartz,  in  (6| 
and  (7),  we  have  included  nonlinear  elastic  terms  only  and 
kept  the  electric  and  electroelastic  terms  linear,  in  accor¬ 
dance  with  the  discussion  in  the  Introduction.  In  Eqs.  (6),  (7), 

and  ( 10)  is  the  electric  scalar  potential  and  Cur.«s>  ^u/icAr4s> 

^LMABNKUf  ^BiMt  *<><1  ^  "c  tiw  second-ordcT  elastic,  third- 

order  elastic,  fourth-order  elastic,  piezoelectric,  and  dielec¬ 
tric  constants,  respectively.  The  associated  boundary  condi¬ 
tions  at  the  surface  of  the  piezoelectric  solid  are  on  or 

u^,bbANlDl  or  g7,  where  IVy,  denotes  the  unit  normal  to  the 
undeformed  surface. 

A  schematic  diagram  of  the  piezoelectric  rod  is  given  in 
Fig.  1,  where  21>2w>5h  and  the  surfaces  normal  to x,  are 
fiiUy  ekctroded.  Since  we  are  interested  in  solutions  which 
vary  slowly  compared  to  2A,  to  lowest  order  the  electrical 
variables  may  be  written  in  the  form*’^ 

=  (11a) 

r,--F/2A,  (lib) 


M 


FIG.  1.  Schematic  diasram  of  the  dectroded  piezoelectric  rod. 

/>,  =  Z),(jr„r),  (lie) 

where  V  is  the  voltage  across  the  electrodes,  and  Z>,  and  D, 
are  not  needed  here,  and  to  this  order  deviations  from  what  is 
ihown  in  Eqs.  (1 1)  are  negligible.  Since  the  electrical  behav¬ 
ior  is  linear  and  the  linear  electric  constitutive  equation  ad¬ 
justed  for  the  rod  is  well  known,"  it  is  convenient  to  omit  the 
electric  field  from  (7)  while  obtaining  the  nonlinear 
rod  equation  and  then  introduce  the  known  electric  quanti¬ 
ties  suitably  adjusted  for  the  rod." 

Since  we  are  interested  in  wavelengths  of  the  order  of  / 
and,  hence,  solutions  varying  slowly  compared  with  A  (and 
w),  both  the  inertia  and  stillhess  of  the  electrodes  are  neg^- 
ble  and  with  the  rod  axis  along  ff,  we  have  the  boundary 
(xmditions 

if;,=0  at  Jrj=  ±w,  ffjy  =  0  at  jr,  =  ±  A.  (12) 

Since  for  the  low  frequencies  of  interest  the  solution  varies 
slowly  compared  to  the  cross-sectional  dimensions  and  flex¬ 
ure  is  unGOUjfled  from  extension,  which  we  are  considenng, 
the  transverse  accelerations  are  negligible  exactly  as  in  the 
linear  theory  of  elastic  rods  and  we  have 

iijsO,  iijsrO,  (13) 

while  we  must  retain  the  displacement  gradients  »2.m  *>id 
Uijt  that  quasistatically  accompany  the  dynamic  longitudi¬ 
nal  deformation  W|,i .  In  view  of  Eq.  (12),  the  aforementioned 
slow  variation  and  the  smallness*  of  w  and  A,  we  see  that  we 
should  take 

Ar,,=0,  ATjy  =  0  everywhere.  (14) 

Furthermore,  from  Eqs.  (1),  (13),  and  (14),  we  have 

ffiz.i  =0,  ff„.,  sO.  (15) 

and  hence  in  the  nonlinear  description  being  obtained  here 
AT, 2  and  K,,  should  also  vanish  approximately  because  they 
vanish  on  the  surfaces  at  AT,  =  ± /.  Sinceyy  ^  is  nonsingular, 
from  Eqs.  (S)  and  (14)  we  have 

Tjm  =0,  Tyff  =0,  (16) 

which  shows  that  in  this  description  (ff  long  wavelength  ex- 
tensional  vibrations  of  elastic  rods  all  Tut,  vanish  except  r,,. 
At  this  stage  the  only  nonzero  components  of  the  Piola- 
Kirchhoff  stress  tensor  are  the  Ar,y,  which  from  Eqs.  (4) 
and  (5)  are  given  by 

A:„  =  (i-h«,,,)r„.  (17) 

^12  ~  *<2,1  T'li,  ff|j  =  W3J  Til.  (18) 

In  the  approximation  being  made  here  u  1 1  is  large  and  all  the 
other  Ul,m>  which  arise  because  of  the  vanishing  of  the 
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{LM  #ll),aFeexpressibleintermsofu,,,  and  are  somewhat 
smaller.  These  relations  are  derived  in  the  next  section  and  in 
the  course  of  the  derivation  it  is  shown  that  for  wavelengths 
large  compared  to  the  cross-sectional  dimensions  of  the  rod 


Mj.,  =0.  «3.,  sO.  (19) 

Thus,  from  Eqs.  (18)  and  (19)  we  have 

/:,2=o,  a:„=j0,  (20) 

and  from  Eqs.  (14)  and  (20)  we  see  that  all  =  0  except 
Jtr,,.  Under  these  circumstances  all  that  remains  of  Eq.  (1)  is 

AT,,.,  =Ai.  (21) 


which  is  the  basic  differential  equation  for  the  anisotropic 
rod.  We  must  now  find  the  nonlinear  relation  between  K^^ 
and  u,.,  that  results  from  Eqs.  (16).  This  is  done  in  the  next 
section. 


The  iterative  procedure  consists  of  first  solving  the  lin¬ 
ear  portion  of  Eq.  (26)  for  „r„  where  the  zero  denotes  the 
order  of  the  iterate,  and  then  substituting  into  the  linear  por¬ 
tion  of  Eq.  (27)  to  obtain  oEr-  Then  the  substituted 
into  the  quadratic  terms  in  Eq.  (26)  to  obtain  ,r„  which  is 
then  substituted  into  Eq.  (27)  along  with  the  substitution  of 
oEj.  into  the  quadratic  terms  to  obtain  Finally,  the  ,£7- 
is  substituted  into  both  the  quadratic  and  cubic  terms  in  Eq. 
(26)  and  only  quadratic  and  cubic  terms  in  £|  are  retained  to 
obtain  jT,.  Accordingly,  the  zeroth  iterate  takes  the  form 

oF)  ~  Ef/Sii,  (28a) 

oEj- =  ST]Ef/Sfi,  (28b) 

which  are  the  well-known  results  in  the  linear  theory.  Now, 
substituting  Eqs.  (28)  into  Eq.  (26)  and  retaining  terms  no 
higher  than  quadratic,  we  obtain 


II.  NONLINEAR  ANISOTROPIC  ROD  EQUATIONS 

Since  all  stress  components  Tlm  have  been  taken  to 
vanish  except  T,  „  Eq.  (7)  can  be  put  in  a  more  useful  form  for 
our  purposes.  To  this  end  we  fiist  write  Eq.  (7)  in  the  usual 
compressed  matrix  notation,  thus 

fp  =  Cf^  Eq  —  if  ic  ^  CpQn  Eq  Ep 

+  i  ^roKS  Eq  E/f  Eg,  (22) 

where"* 

P.Q.R,S,T=  1,2 . 6.  (23) 

and  we  recall  from  Eq.  (9)  that  the  Eq  are  nonlinear  in  the 
Vjr.  L  ■  now  operate  on  Eq.  (22)  with  the  reciprocal  matrix 
of  the  second-order  elastic  constants  cf, ' ,  which  ate  the  elas¬ 
tic  complianoes  Sj-p,  to  obtain 

Et  *  •*ri  4-  dgp  ^ g  —  ^  ijp  CpQg  Eq  Eg 

—  I  i  ff  Eq  Eg  Eg,  (24) 

where  we  have  employed  the  condition  that  all  vanish 
except  r,  and 

^KT  *  (25) 

is  a  defined  piezoelectri  j  constant  that  is  particularly  useful 
when  approximatioiu  ate  made  as  a  result  of  traction-free 
boundary  conditions.  For  purposes  of  this  work  it  is  conven¬ 
ient  to  write  the  six  equations  in  Eq.  (24)  as  two  separate 
equations,  one  for  r=  1  and  the  other  five  for  r  ^  1,  which 
in  the  absence  of  the  piezoelectric  term  dgr  9 g  take  the 
forms 


|7’|  —  Ey/s^  I  -f  (l/2rJi  )ri/>  <^pQg  s^Q  s^g  E  J ,  (29) 

which  is  the  first  iterate  for  T,  and  is  all  we  would  need  if  we 
want  a  quadratic  relation  for  .  Substituting  from  Eq.  (29) 
and  Eq.  (28b)  into  Eq.  (27)  and  retaining  terms  no  higher  than 
quadratic  in  £,,  we  obtain 


\Et  J  ~  ^Q\  ^R\  E\, 

r,,  2st,  \s,, 

7*  9^1-  (30) 

The  further  substitution  of  Eq.  (30)  into  Eq.  (26)  and  the 
retention  of  terms  no  higher  than  cubic  in  £|  yields 

^T^=cE^-\-E^E\  ■^Y\E\,  (31) 

where 


^R\> 

■*11  ^  ■^1  ^ 

1  [  Aei  ^ 


+  ~ i\p  CpQgg  jg ,  Sgi  is  1 J,  (32b) 

where  the  caret  on  the  subscript  Q  means  that 

jgK=0fore=l.  (33) 


The  second  iterate  for  T,  given  in  Eq.  (3 1 )  with  Eqs.  (32)  is  the 
cubic  relation  we  need  to  obtain  the  cubic  constitutive  rela¬ 
tion  for  £, ,  in  terms  of  u,,, .  Substituting  from  Eq.  (31)  with 
Eq.  (9),  and  Eq.  ( 1 9)  into  ^.(17)  and  reintroducing  the  linear 
piezoelectric  term,  we  obtain  the  result 

A’m=c«,,,  -f-£(«,.,)*-f-y|i/,.,)^-(</3,A,,)8p3,  (34) 


where 


T,=  -LE,+lflLc  EgEg 

■*11  2  5,1  J 

"!■  ~z  ^PQRs  Eg  Eg  Eg,  (26) 

6  5,,  * 

Et  —  ~  i  ^TP  <^pQR  Eq  Eg 

~  l^TP  (pQRs  Eq  Eg  Eg,  T^\.  (27) 


E  =  +  (35a) 

r  =  ri  +  y3,  +  ic.  (35b) 

Equation  (34),  with  Eqs.  (35)  and  (32),  is  the  constitutive  rela¬ 
tion  we  have  been  seeking  for  the  arbitrarily  anisotropic  pie¬ 
zoelectric  rod,  which  is  nonlinear  in  the  mechanical  dis¬ 
placement  gradient  u,,,  but  linear  in  the  electric  field. 

We  must  now  show  that  £,2  and  £„  do  indeed  vanish 
approximately  as  indicated  in  Sec.  I.  To  this  end  we  note  that 
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(47) 


at  the  zeroth  iterate,  from  Eq.  (28b)  with  Eq.  (9),  in  tensor 
notation,  we  have 

='•*1.11  “i.i.  (36) 

where 

Skl  =  +  “/,*)»  (37«) 

'Viii  (37b) 

and  u,  =  u,(Jir„  r).  For  KL  =  22  and  33,  respectively,  we 
have 

“jj  = '’2211  “l.l*  *<3.3  =  ^3311  “1.1  •  (38) 

which  may  be  integrated  immediately  to  give 

•*2 —•^2^2211  **1.1  ■l"/'(^l>'^3>  *)> 

*<3  “■3^3^3311  **1.1  +  8(-^l» -3^2*  ^  )•  (39) 

From  the  remaining  three  relations  in  Eq.  (36)  with  Eq.  (37a), 
we  have 

i(**l.2  +  **2,1 )  ==  *’1211  **l.l  >  (40®) 

J(**i.3  +  **3.1 )  ~  *‘1311  **1.1 »  (40b) 

J(**2.3  +**3.2)  =  '’2311  **1.1-  (40c) 

Since  the  rod  axis  lies  along  and  we  have  assumed  that 
torsion  is  uncoupled  from  extension  in  this  thin  anisotropic 
rod,  we  have 

**'32  =  i(**2.3  -**3.2)  =  0,  (41) 

where  W32  is  the  small  (linear)  local  mechanical  rotation 
about  the  rod  axis  AT,.  From  l^s.  (39),  (40c),  and  (41),  we 
obtain 

/= -^3^2311  1*1,1  •  J?  =  •^2*'2311  **1.1  •  (42) 

Now,  substituting  from  Eqs.  (39)  and  (42)  into  (40a)  and 
(40b),  respectively,  we  obtain 

**1.2  =2r, 211 1*1,1  —  (Ar2r22n +A^3*’23ii)**i.ii, 

(43) 

**1,3  =  2^1311  **1.1  ~  {^3^3311  ■|■•^2*■231l)**l.ll■ 

Since  we  are  interested  in  the  limit  of  long  wavelengths  (A )  as 
the  rod  collapses  to  a  line,  we  have 

maxIATal  lu  ^  =  (44) 

A  A  Ai  Ai 

which  with  Eqs.  (39),  (40c),  (41),  and  (43)  yields 


o**i,2  —  2ri2ii  *<|,|,  (45a) 

o**M  =  2ri3ii  Mi  l,  (45b) 

o**2.i  =0,  (45c) 

o**2,3  ='■2311  **1.1.  (45d) 

o**3,i  =0,  (45e) 

o**3a  =  *’2311  **1.1 »  (43f) 


which  along  with  (38)  gives  all  the  zero-order  displacement 
gradients.  Equations  (45c)  and  (45e)  show  that  at  the  zeroth 
iterate  Eqs.  (19)  hold  and  since  ,Ktj  is  obtained  from  o**.4.«> 
we  have  shown  that 

lATijsO,  ,Ar,3=0.  (46) 

In  the  case  of  the  first  iterate,  from  Eq.  (30)  with  Eqs.  (9) 
and  (45),  in  place  of  Eq.  (36),  we  obtain 
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where  in  the  compressed  matrix  notation 


Xtfoii  Jgi  3*1 -l-<*ri»  T ^l,  (48) 

and  each  On  consists  of  a  sum  of  terms,  each  of  which  con¬ 
tains  a  quadratic  combination  of  the  rn  from  Eqs.  (45)  and 
(9).  The  exact  forms  of  the  a-n  are  cumbersome  to  write 
because  of  the  forms  in  Eqs.  (45)  and  are  not  needed  here. 
From  Eqs.  (47)  and  (37b)  along  with  the  procedure  followed 
in  obtaining  Eqs.  (45)  in  the  case  of  the  zeroth  iterate,  it  is 
clear  that  we  can  obtain  relations  similar  to  those  in  (45)  if  we 
can  obtain  the  relation  equivalent  to  Eq.  (4 1 )  in  the  case  of  the 
first  iterate.  It  is  shown  in  the  Appendix  that  to  second  order 
in  the  mechanical  displacement  gradients  the  local  or¬ 
thogonal  transformation  tensor  is  given  by 

^KL  —^KL  (49) 

where 

Hlk  =  **'ijt  +  i  U>LM  ~  ^LM  ^Mk)> 

(50) 

and 

^uc  ~  J(**it  j.  ~  **t,x)'  (31) 

In  the  linear  theory  the  same  type  of  procedure  shows  that 

Rkl—S/cl-^Wu^.  (52) 

Since  in  the  quadratic  case  iJ/jc  takes  the  place  of  W/jc  in  the 
linear,  1123  represents  the  torsional  rotation  of  the  rod  in  the 
quadratic  case.  Consequently,  in  order  to  uncouple  torsion 
we  take 

1223  =  0.  (53) 

Substituting  fiom  the  zeroth  order  results  in  the  nonlinear 
terms  in  Eq.  (50)  and  using  the  limits  in  Eqs.  (44),  from  Eq. 
(53)  we  obtain  for  the  first  iterate 

•*^23  =  i  *’l2ll*'l3ll(**l,l)^  (54) 

in  place  of  Eq.  (41)  for  the  zeroth  iterate.  From  the  forms  in 
Eq.  (47)  and  Eq.  (54),  which  take  the  place  of  Eqs.  (36)  and 
(41),  respectively,  it  is  clear  that  by  foUowing  the  procedure 
used  earlier  in  the  case  of  the  zeroth  iterate,  we  obtain  ex¬ 
pressions  similar  to  those  obtained  with  the  zeroth  iterate 
but  containing  additional  terms  depending  on  (u,,,)^.  Fur¬ 
thermore,  when  we  take  the  limit  in  Eqs.  (44)  we  obtain 

i**2.i  =  0,  iM3,|  =  0,  (55) 

and  since  2^^  is  obtained  from  from  Eq.  (18)  we  have 
shown  that 

2^12~0,  2^|3S;0.  (56) 

When  all  stress  components  Tijc(K^)  vanish  except 
T\  i(Ar, ,)  the  linear  piezoelectric  constitutive  equation  can  be 
written  in  a  more  useful  form.  To  this  end  we  first  write  Eq. 
(6)  in  the  compressed  matrix  notation,'*’  thus 

2)z.=*^L0^0+«ix  (57) 

and  we  substitute  from  the  linear  portion  of  Eq.  (24)  into  Eq. 
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from  which,  with  the  aid  of  Eqs.  (11),  we  obtain 

^3=-^«u+(eri-— Vi.  (61) 

•*1111  '  ■*1111' 

in  which  we  have  reintroduced  the  tensor  notation. 

The  substitution  of  Eq.  (34)  with  Eq.  (lib)  into  Eq.  (21) 
yields 

cai.ii  +/9I(«i.iPl.i  +r[(«i.i)*l.i  =/>«i.  (62) 

which  is  the  differential  equation  of  motion  of  the  piezoelec¬ 
tric  (anisotropic)  rod.  The  boundary  conditions  consist  of  the 
specification  of  ATn  or  u,  at  the  ends  of  the  rod  JIT,  =  ±1. 
After  a  solution  has  been  obtained,  the  current  /can  be  cal¬ 
culated  from  the  relation 


'  =  2uij^  /)3 


and  we  recall  that  in  Eqs.  (34)  and  (61)  is  related  to  the 
driving  voltage  Kby  Eq.  (11b). 

The  quadratic  rod  coefficient  p  has  been  calculated 
from  Eq.  (3Sa)  with  Eqs.  (32a)  for  a  number  of  orientations  of 
quartz  rods  using  the  second-"  and  third-‘^  order  elastic 
constants  of  quartz  and  the  results  are  plotted  in  Figs.  2-4. 
The  linear  rod  constant  c  is  plotted  ak^  with  — ^ineach 
figure.  The  quadratic  rod  codficient  P  is  always  negative  and 
ranges  in  value  from  about  —  ctoabout  —  Sc.  Figure  2  is  for 
the  rotated  X  cuts  with  $  =  V  corresponding  to  the  rod  axis 


no.  3.  Linear  and  quadratic  elastic  constants  for  rotated  K-cut  quartz 
rods. 

along  the  crystallographic  Y  direction.  The  locations  of  the 
-H  5*  — ,  hfT  —  ,  and  —  18*  —  JT  cuts  are  shown  in  the  fig¬ 
ure.  Figure  3  is  for  the  rotated  Y  cuts  with  0  =  (f  corre¬ 
sponding  to  the  rod  axis  along  the  crystallographic  X  direc¬ 
tion.  This  figure  exhibits  the  twofold  symmetry  of  quartz 
with  X  the  digonal  axis.  Figure  4  is  for  doubly  rotated  cuts 
with  6  =  V  corresponding  to  the  rod  in  the  X-Y  plane  at 
-t-  4S*  from  the  axis.  In  Figs.  2-4, 6=  ±  90*  corresponds 
to  the  rod  axis  along  the  crystallographic  Z  direction. 


III.  INTERMODULATION  IN  RODS 

In  this  section  we  consider  the  problem  of  intermodula¬ 
tion  in  piezoelectric  rods.  For  clarity  we  reproduce  certain  of 
the  pertinent  equations  obtained  in  Secs.  I  and  II  here.  The 


FLO.  2.  Linear  and  quadnuk  elaitk  conttanto  for  routed  2r<ul  quartz 


FLO.  4.  Linear  and  quadratk  elastic  constants  for  a  family  of  doubly  routed 
quartz  rods. 
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stress  eqiwticMi  of  motion  and  charge  equation  of  electrostat¬ 
ics  for  the  piezoelectric  rod  take  the  form 

A’n.,  (64a) 

!),,,=  0,  (64b) 

where 

=  c«,.,  -h/ff  («,.,  f  +  y(«,,,  ^3.  (65) 


/)j  =  (rf3,/s„)tt,.,  -|-(f3^j  -dlt/Sn)^i,  (66) 

and  c,  P,  and  y  arc  given  in  Eqs.  (32)  and  (35).  From  Eq.  (11b) 
we  have 

gp3=  -(F/2A )«'“',  (67) 

where  V  is  the  magnitude  of  the  applied  voltage  and  a  is  the 
steady  driving  frequency.  For  a  rod  free  at  both  ends  the 
boundary  conditions  are 

A’,,  =  0  at  AT,  =  ±/.  (68) 

Since  the  asymptotic  solution  to  the  nonlinear  problem 
is  obtained  by  iteration  from  the  solution  of  the  appropriate 
linear  problem,  we  begin  the  analysis  with  a  brief  presenta¬ 
tion  of  the  known  solution  of  the  associated  linear  problem. 
As  a  solution  of  the  linear  problem,  i.e.,  Eqs.  (64a)  and  (68), 
with  Eq.  (65)  in  the  absence  of P  and  y,  we  write 

oii,  =  lott-(</„/2A)KY,y“*',  (69) 

the  substitution  of  which  in  Eqs.  (64a)  and  (68)  yields 

Co«.„  =/>V(d3i/2A  )FAr„  (70) 

ott.,  =0  atA',=  ±/.  (71) 

The  eigensolutions  of  the  homogeneous  form  of  Eqs.  (70)  and 
(71),  i.e.,  with  K  =  0,  may  be  written  in  the  form 

Ui„=u„e  (72a) 

u„  =a„  sin  ri„X,.  (72b) 

Substituting  from  Eqs.  (72)  into  the  homogeneous  form  of 
Eqs.  (70)  and  (71),  we  obtain 

p°al=crfl,  (73a) 

Ti„  =  nir/2l,  n  odd.  (73b) 

We  write  the  steady-state  solution  of  the  forced  vibration 
problem  in  the  form 

o»  =  X^"  S’"’?"'*’" 

It 

the  substitution  of  which  in  Eq.  (70)  along  with  the  use  of  the 
orthogonality  of  Eq.  (72b)  with  Eq.  (71)  and  Eq.  (73a)  yields 

A.  =  -  8/(  -  1)'-  '  V/\  [{col/w^)  -  1  ]nV2h  1 . 

(75) 

Since  Eq.  (75)  has  a  resonance  denominator,  in  the  vicinity  of 
a  resonance,  say  the  TVth,  one  term  in  the  series  in  Eq.  (74) 
dominates  and  with  Eq.  (69)  we  have 

0«,  =  [  sin  {NirX,/2l )  -  (dj,  /2A  )X,  ] (76) 

where  we  have  introduced  the  notation  <u„  and  ,Af^  to 
denote  the  quantities  at  driving  frequency  a),.  Substituting 
from  Eqs.  (76)  and  (67)  at  <u,  into  Eq.  (66),  which  is  then 
substituted  into  Eq.  (63),  we  obtain 


L,  =  Awlic,  1  -h  -  ,  (77) 

where  we  have  made  the  identification  =  —  I  ^X(o^  since 

at  <U|  the  crystal  is  passive'^  and 

(^3l)*  (78) 

and  we  have  replaced  as  by  an  where 

(79) 


in  which  is  the  unloaded  quality  factor  of  the  resonator  in 
that  mode,  and  serves  to  prevent  the  resonant  denominator 
from  vanishing  ^1(0  =  0}^. 

Since  we  are  interested  in  intermodulation  due  to  the 
simultaneous  application  of  test  tones  to  the  resonator  at  two 
nearby  frequencies  a ,  and  <U2,  both  of  which  are  in  the  neigh¬ 
borhood  of  &>jv,  we  consider  the  steady-state  linear  solution 
to  the  vibration  problem  already  presented  in  this  section  to 
be  at  the  frequencies  a ,  and  Moreover,  since  we  are  inter¬ 
ested  only  in  the  situation  in  which  the  intermodulation  fre¬ 
quency  il  is  in  the  vicinity  of  only  certain  combinations 
of  the  nonlinear  frequency  products  are  significant.  To  ob¬ 
tain  expressions  for  the  nonlinear  products  in  Eq.  (65)  we 
first  write  the  dominant  real  solution  in  complex  notation, 
thus 


«1  =i(o«>-I^“''  +  0«Ar2«'“'' 

+  +  (80) 

where  *  denotes  complex  conjugate  and  we  have  ignored  the 
terms  containing  and  because  they  do  not  have  reso¬ 
nant  denominators, and  then  substitute  into  the  quadratic 
product  in  Eq.  (65)  to  obtain  a  linear  steady-state  inhomo¬ 
geneous  forced  vibration  problem  for  the  first  iterate  solu¬ 
tion.  When  the  appropriate  linear  inhomogeneous  vibration 
problem  has  been  defined  standard  complex  notation  for  lin¬ 
ear  systems  is  employed.  This  type  of  procedure,  i.e.,  the  use 
of  the  previous  iterate  to  define  the  linear  inhomogeneous 
problem  for  the  next  iterate  in  the  above-mentioned  manner, 
is  employed  at  each  successive  stage  of  iteration.  In  view  of 
the  foregoing  statements,  the  form  of  the  linear  solution,  and 
previous  experience'^  only  the  frequencies  of  (2<u,  -  ruj)  and 
(2<t>2  —  <»i)  are  of  interest  in  our  intermodulation  study.  Fur¬ 
thermore,  it  is  sufficient  for  our  purposes  to  determine  the 
response  at  (2<u,  — Oj)  only  because  the  response  at 
(2a>2  ~  )  can  then  be  obtained  simply  by  an  interchange  of 

the  subscripts  1  and  2. 

Since  we  are  interested  in  the  response  at  (2<u,  —  (Uj) 
only,  the  form  of  Eq.  (65)  along  with  (80)  and  the  sur¬ 
rounding  discussion  reveals  that  we  need  determine  the  first 
iterate  solution  due  to  the  quadratic  nonlinearity  only  at 
2<U|  and  (tu,  —  Furthermore,  in  the  interest  of  brevity  we 

obtain  ,u,  at  2<ii>,  only  and  simply  present  the  changes  in  the 
result  at  (02 ,  —  As  already  noted  the  first  iterate  equation 
at  2<U|  is  obtained  by  substituting  gUt  itt  (76)  in  the  qua¬ 
dratic  term  in  Eq.  (64a)  with  Eq.  (65),  in  accordance  with  Eq. 
(80)  and  the  associated  discussion,  but  ignoring  and  the 

time-independent  term  and  neglecting  the  term  without  a 
resonant  denominator'^  to  obtain 


=^,Ai^j7^(sin2?7^Ar,/2)e''“'’.  (81) 


2007 


J.  Acoutt.  Soc.  Am.,  Vol.  73,  No.  6,  June  1963 


H.  F.  Tiersten  and  A.  Ballato:  Vibrations  of  rods 


2027 


Since  the  boundary  condition  is  given  by  Eq.  (68),  by  virtue  of 
Eqs.  (73)  the  solution  for  ,u,  takes  the  form 

,«,  =  (X  sin  Jr,  -  AT, V”"’'.  (82) 

\  m 


where  the  B„  are  obtained  from  the  orthogonality  of  the 
sin  i7„  Xy  Including  only  one  term,'^  say  the  ^th,  in  the 
series  in  Eq.  (82),  we  have 

“  2*i“(4s)i -3i,|  [(Si,/*!ii|-ll«V2*’ 


where  we  have  used  Eq.  (73a)  and 

_  2//'sin(27V-itf)ff/2  sinjlN  +  M)ir/2\  .... 

"A  2N-M  2N  +  M  )' 

Substituting  from  Eqs.  (82)  and  (67)  at  2<u,  into  Eq.  (66), 
which  is  then  substituted  into  Eq.  (63),  we  obtain 

r  _  ^33/,  , 


IiLu.  =  —  ^li2o> 


2h\ 


- ■2AV  ' 

5„ 


ili! _ V 


where  we  have  made  the  identification  =  /  at  2<u,  since 
at  2n>,  the  crystal  is  active.'^  In  a  similar  way,  again  includ¬ 
ing  only  the  A/th  term,'^  at  (<u,  —  0)2)  =  a,  we  obtain 

,ii,  =  (c,,  sin  r,^  X,  -  jr.V"'.  (86) 


P\  *Vn  PnM  _  8/(—  1)‘*^  **'^<^31^0 

VfP[oP^  —  Sv )  [  (SJ,/a^)  —  1  ]  ’ 


I„  =  -4w/io^^l  + 


8(W__V 

(o>l,/a^)  -  1  ]  / 


2h\ 

(88) 

s„  /pV’-SJ,) 

The  linear  second  iterate  equation  at  /2  =  (26),  —  6)2)  is 
obtained  by  substituting  the  appropriate  expression  for  ,i<„ 
i.e.,  a  linear  combination  of  the  zeroth  and  first  iterate  solu¬ 
tion  functions  at  6>„  6);,  26>„  and  a,  which  is  of  the  form 

i«i=  -  -~X,Je’“'' 

+  sin  Vjy  - 

,  ( n  ■  V  <^31^2®,  \  fl®,, 

+  sm  ijin  X I  —  — — —  X I  je 

+  (c^  sin  X,  -  +  c.c.],  (89) 

in  both  the  quadratic  and  cubic  terms  in  Eq.  (64a)  with  Eq. 
(6S),  neglecting  terms  without  resonance  denominators,  i.e., 
terms  containing  F,  and  I2  hut  not  Pj®,  nnd  compared 
to  terms  with  resonance  denominators  and  retaining  all 
terms  of  order  ,A  j,2^  %  and  proportional  to  e"*',  with  the 


f2«l.ll  -P  2«1 

=  (^/2)[(2/f  *  \A,^Cm)ffN  Vm 
X(i7Arsini7>r  X,cosrjM  X, 

+  Vm  sin  Vm  cos  Ti,^  Xt) 

~  ’7y(^3i/2A  )(2/4  *  1^2®,  +  [Ay  )sin  V.v 
+  S  ViA  ItA  I  ri%  cos^  i}y  Xi  sin  27.V  X,e‘^'.  (901 

As  in  the  case  of  the  zeroth  iterate  we  take  a  series  solution 
for  jU,  in  the  form 

2«,  =  (X  sin  r,„X,  -^Vn  (91 ) 

and  with  the  aid  of  the  orthogonality  of  the  sin  't]„Xy  we  find 
for  n  =  N,  which  dominates  the  others 

^  ^  J_ 

""  [(S^//2")- 1]JW2A'^  llp'^ 


^iiA  *  B^  -(-  iAfj  Cu)'!). 


N  Vm(jJn 


Vn9nm  + 


Vm  Pnm 


Pd2\  Vn^^A  %  ^2®,  +  tA\  ) 
4h  p°({}^—a%) 

9  YiA  \  *Vn 
le  p<^in^~^i)’ 


where 


VsM  — 


sin  jfJ 


_  1  ^sin(2);jv 

2V  2vy- 


+  Vm)I  ,  M^Vn-Vm)! 


'n  -  VAr)/\ 
,  -  Vm  ) 


2\  2i7a-  +  Vm  2)7^  -  Vm  ) 

Substituting  from  Eqs.  (91)  and  (67)  at  fl  into  Eq.  (66),  which 
is  then  substituted  into  Eq.  (63),  we  obtain 

r  A...un^n(y  .  8(^  )^  N 


2h  \  7r‘'N^\(o>%/n  ^)  —  1  ]/ 

j,i  p'‘co%[(n yo>%)  - 1]  V 


j,i  p°co]j[(n^/o}]^)  - 1]  V 
+  xAyC^]VsVM^  -  ^^(2^  *  +  ^As 


—  in  2w  —  —  ^  iAi2A%V\(— 

5,1  16  p°a}%  [{n~/a%)—l] 


where  we  have  made  the  identification  1^,=!  dxn  since  the 
crystal  is  active'^  at  O  and 

fNM  =  VnVnm  +  Vm  Pnm/2.  (95) 

A  schematic  diagram  of  the  circuit,  which  is  driven  at 
the  two-test  tone  frequencies  o),  and  6)3,  is  shown  in  Fig.  5, 
where  y^  is  the  driving  voltage,  Rg  the  generator  resistance, 
and  /?2,  the  load  resistance.  Application  of  KirchhofTs  vol¬ 
tage  equation  to  the  circuit  shown  in  Fig.  5  at  6),  and  6)1 
yields  the  two  equations 

.Vg+I^^(Rg+R,]+y,,^=.0, 

2y,+/.JBg+R,)+y^^=0, 

which  with  Eq.  (77)  for  6),  and  the  equivalent  for  6),  enables 
the  determination  of  and  in  terms  of ,  F  and  jFL, 
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(100) 


FIO.  S.  Schematic  diagram  of  the  reduced  test  circuit. 


respectively,  from  which  and  act  obtained  in  terms 
of ,  F,  and  2F,,  respectively,  from  Eq.  (75).  The  application 
of  Kirchhoff  s  voltage  equation  to  the  circuit  shown  in  Fig.  3 
at  2tu„  a,  and  Q  yields 

+  ^2.)  +  >;  =  0,  =  2ta„  a,  /3,  (97) 

since  at  th^  frequencies  K,  =  0.  From  Eq.  (97)  with 
K  =  2ca,  and  K  =  a.  respectively,  and  Eq.  (83)  and  Eq.  (88), 
Vjt,,  F„  may  readily  be  determined  in  terms  of ,  F,  and 
jF,.  Finally,  from  Eq.  (97)  with  k  =  /J  and  Eq.  (94),  F^  can 
readily  be  determined  in  terms  of  |  F,  and  2Ff  if  is  known. 
However,  since  actually  x  is  the  only  unknown  quantity,  a 
measurement  of  the  load  voltage  enables  the  deter¬ 
mination  of  la  nnd  Va  from  Eqs.  (97),  from  which  y  csn  be 
evaluated  frtm  Eq.  (^)  for  any  desired  orientation,  and  in 
particular,  for  the  practical  3*  X  cut. 


IV.  NONUNEAR  RESONANCE  IN  RODS 


In  this  section  we  consider  the  problem  of  nonlinear 
resonance  in  piezoelectric  rods.  The  basic  nonlinear  equa¬ 
tions  required  for  the  problem  of  interest  here  are  Eqs.  (64)- 

(68) .  We  again  obtain  the  solution  by  iteration  from  the  linear 
solution,  but  now  the  successive  iterates  must  be  corrected 
because  we  are  interested  in  the  steady-state  nonlinear  solu¬ 
tion  at  the  driving  frequency  a,  in  which  there  is  a  nonlinear 
dependence  of  the  resonant  frequency  on  the  amplitude 
of  the  extensional  displacement  of  the  r^.  Clearly,  the  ze¬ 
roth  iterate,  or  linear  solution,  is  given  in  Sec.  Ill  and  Eqs. 

(69) -(79)  hold  for  purposes  of  this  section.  Again,  in  the  vi¬ 
cinity  of  the  fVth  restmance  one  term  dominates  the  series 
and  the  zeroth  iterate  is  given  by  Eq.  (76)  with  the  subscript  1 
stricken  from  thei^jy  and  the<u. 

In  rccordance  with  the  discussion  surrountting  Eq.  (80) 
in  order  to  obtain  the  linear  equation  for  the  first  iterate,  the 
zeroth  iterate  solution  is  written  in  the  form 


[(i^jvsin^jr,- 


(98) 


which  b  substituted  in  the  quadratic  term  in  Eq.  (64a)  with 
Eq.  (63)  to  obtain  the  linear  inhomogeneous  equation  for  the 
first  iterate 


n  sin  217^^  XM  j,^  +  Af,A  J). 

(99) 


Since  the  boundary  condition  b  given  by  Eq.  (68)  with  Eq. 
(63)  and  the  Fi  b  negligible,  by  virtue  of  Eq.  (73b)  the  solu¬ 
tion  for  ,11,  takes  the  form 


C„  sin  Vm  Xt. 

m 

where  the  B„  and  C„  are  obtained  from  the  orthogonality  of 
the  sin  ri„  Xi  and  we  have  made  use  of  the  fact  that  the 
electrodes  are  shorted  at  zero  frequency.  Since  the  case  of 
primary  interest  is  TV  =  1  and  for  that  case  M  =  1  gives  the 
largest  term  in  each  series  in  Eq.  (100)  and  also  results  in  the 
largest  single  term  caused  by  the  series  in  Eq.  (100)  at  the 
next  iterate  and  we  are  not  interested  here  in  an  exhaustive 
treatment  of  all  possibilities,  but  only  in  indicating  the  na¬ 
ture  of  the  influence  of  the  terms  in  the  series  in  Eq.  (100)  on 
the  nonlinear  resonance  equation  and  obtaining  the  expres¬ 
sion  due  to  the  term  of  greatest  influence  for  the  case  =  1 , 
we  consider  the  term  Af  =  TV  only  in  each  series  in  Eq.  (100). 
Accordingly,  including  only  the  A/th  term  in  both  series  in 
Eq.  (100),  with  the  aid  of  the  orthogonality  of  the  sin  AT, 
and  Eq.  (73a),  we  have 


%Vn  Pnm 
2lp'^{4o^  -  0)1,) 

_  nVnPnm 
2/po<s>L 


(-l)l*/-.|/28/<f3,K,„ 
\(o)\,/Ao)^)  —  1]A/ 

(101) 


where  given  in  Eq.  (84).  Substituting  from  Eq.  ( 100) 
and  Eq.  (67)  at  hi)  into  Eq.  (66),  which  is  then  substituted  into 
Eq.  (63),  we  obtain 


-  ^u)li2o)  ^ 
2h 


\[(i)\,/Ao)^)  —  l]TVf  V 


+  i2o,2w^ 


J,,  lp'^(Ao)^  -  0)1,) 


(102) 


where  we  have  made  the  identification  =  /  at  2o),  since 
at  2o)  the  crystal  is  active.'^  In  accordance  with  earlier  db- 
cussions  in  order  to  obtain  the  linear  equation  for  the  second 
iterate,  the  first  iterate  solution  is  written  in  the  form 


1  [/.  .  Nit  „  \ 

+  sin  ^  AT,  - 
+  sin-^AT,  -|-c.c.j. 


(103) 


and  substituted  in  both  the  quadratic  and  cubic  terms  in  Eqs. 
(64a)  and  (68),  with  Eq.  (65),  while  retaining  all  terms  of  order 
A],A^  and  containing  the  time  dependence  e*"'  only  to  ob¬ 
tain  the  homogeneous  differential  equation 


Tj^  AT,  cos  17*,  AT, 

-|-?7*,  sin  t}m  A',  cos?7;v  AT,) 

-  Vs^  % 

-\-^yA],A  X,  sin  17^  Ar,e"“',  (104) 

4 
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and  the  inhomogeneous  boundary  conditions 

s,,  2h 

=  -  y *  Vs  COST! f,  I 

x(5m  VmCOSVmI-  ^  ^20.) 

+  i7,v  cos  ri^  ICj^Vm  cos  17*,  /  je“' 

-  Y  * ’/i  cos^i/jv  at  -y,  =  ±  /, 

(105) 

where  we  have  neglected  terms  without  resonant  denomina¬ 
tors,  i.e.,  those  containing  V^,  compared  to  terms  with  reso¬ 
nant  denominators  and  we  have  not  yet  set  M  =  It  should 

be  noted  that  the  right-hand  side  of  Eq.  (lOS)  turns  out  to  be 
negligible  because  is  small  and  nearly  given  by  Eq. 
(73b).  Since  we  regard  the/^jy,  17/,,,  and  ri^  on  the  right-hand 
sideof  Eqs.  (104)  and  (lOS)  as  presently  unknown  and  ignore 
Eqs.  (73b)  and  (75),  the  equations  are  homogeneous  for 
=  0  and  we  say  the  iterate  is  corrected.  ” 

Before  we  obtain  the  nonlinear  eigensolutions  of  Eqs. 
(104)  and  (105)  we  first  note  that  for  the  short  circuit  case 
^'2.=0,  and  the  open  circuit  case  12^=0  yields  the  greatest 
difference  in  eigenfrequency  from  the  shorted  case.  Further¬ 
more,  the  open  circuit  voltage  obtained  from  Eqs.  (102) 
when  substituted  in  Eq.  ( 104)  contains  the  product  of  (k  j,  )^ 
and  A  %A  %,  both  of  which  are  small.  Since  the  product  of 
two  small  quantities  is  negligible,  we  may  omit  in  Eqs. 
(104)  and  (105)  in  the  sequel.  The  nonlinear  eigensolutions 
must  satisfy  Eq.  ( 104),  wUch  for  M  =  N  takes  the  form 


C2«i.ii  -p5«i  =  sin 2i7Af Xt 


Y'*  s^  sV%  (sin  Vs  X, 


sin  3vs 


where 


H=  -  liLyAiAX^Vs, 

3  yc  I 

J_*L - j,  (107) 

and  we  have  made  use  of  the  fact  that  we  are  obtaining  the 
steady-state  solution  for  ea  in  the  vicinity  of  <ay  only  and  we 
have  employed  Eqs.  (73),  (84),  and  (101).  The  homogeneous 
boundary  conditions  are  given  by  Eq.  ( 105)  with  =  0.  As 
the  iVth  eigensolution  of  the  linear  differential  equation 
(106),  we  take 

,H,  =  (108) 

where 

Uy  =  Ay  sin  Vs  +  ^s  sin  Ivs  -^i  +  ^s  sin  3vy 

(109) 

the  substitution  of  which  in  Eq.  (106)  yields 

=  cv^y(\  +  H),  (110) 


Gs=^K(fs/Nir)pAy, 
Ly  =  —  pAy/%, 


(111a) 

(111b) 


where 


K  =  pycY,  (112a) 

~  ^)V%  ^s^  *  •  (112b) 

Since  the  amplitude  Ay  is  very  small,  Eqs.  (Ill)  with  Eq. 
(112b)  show  that  Gy  and  Ly  are  negligible  for  our  purposes 
except  insofar  that  they  might  result  in  a  change  in  Ss  by 
virtue  of  the  boundary  conditions  (105).  Substituting  from 
Eqs.  (108)  and  (109)  into  Eq.  (105),  we  obtain 

/  (A  fs 

cosvsl  =  -p\cosVsl+  ^ 

V  81  Nir 

- -\-OOs2Vy  1) 

n 


- ^cos  3j7jv  /Y 

24  / 


The  roots  of  Eq.  (113),  with  Eq.  (110),  determine  the  ampli¬ 
tude  dependent  eigenfrequencies  ay  of  this  piezoelectric  rod 
in  the  nonlinear  case.  In  the  absence  of  the  nonlinear  term/r, 
Eq.  ( 1 1 3)  reduces  to  Eq.  (73b).  Since/r  is  small,  the  roots  Vs  ^ 
may  be  obtained  iteratively  from  the  roots  in  the  linear  case, 
which  are  given  in  Eq.  (73b).  Consequently,  in  the  nonlinear 
case  the  roots  must  differ  from  Nir/1  by  small  quantities,  say 
Ay,  and  we  may  write 

Vy  /  =  Nn/l  Ay,  N  oAA.  (114) 

Substituting  from  Eq.  (114)  into  Eq.  (113),  expanding  and 
retaining  terms  linear  in  Ay,  we  obtain 

=  -  (  -  1)'"'-  '^^s{fy/Nn)p,  (115) 

which  shows  that  as  a  consequence  of  the  existence  of  the 
coefficient  P  of  the  quadratic  term,  the  nonlinear  boundary 
conditions  have  a  small  influence"  on  Vs-  Since <,  1,  from 
Eq.  (110),  with  the  aid  of  Eqs.  (114)  and  (115),  we  have 


Nv(  c  ,  p 

'’"'wj  ('+T 


-l-il" 


128  Kfy 


kfs  \ 


which  gives  the  resonant  frequency  of  the  ATth  mode  of  ex- 
tensional  vibration  of  the  thin  piezoelectric  rod  with  small 
piezoelectric  coupling,  including  the  dependence  on  the  am¬ 
plitude  of  vibration." 

The  solution  to  the  nonlinear  forced  vibration  problem 
at  a  driving  frequency  a  in  the  vicinity  of  ay  may  now  be 
obtained  by  means  of  an  expansion  in  the  eigensolutions, 
while  retaining  the  nonlinear  correction  in  the  dominant  N  th 
eigensolution  only.  Accordingly,  we  write 

,«,  =  [,«-(d„F„/2A)Ar,]e‘“',  (117) 

where 

2«  =  Ay  sin  Vs  +  Gy  sin  2vy  X,  -|-  Ly  sin  3vy  X, 

A„  sinVn  X„  «  odd,  (118) 

and  Gy  and  Ly  are  determined  from  Eqs.  ( 1 1 1 ).  In  the  itera¬ 
tive  procedure  we  are  employing  the  influence  of  the  large 
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tenn  in  the  expansion  in  Eq.  (1 18),  i.e.,  ,4^^  sin  -^i.  already  ing  that  resulting  from and  the  first  iterate  solution.  Substi- 

appears  on  the  right-hand  side  of  Eqs.  ( 104)  and  ( 105)  includ-  tuting  from  Eq.  ( 1 1 7)  with  Eq.  ( 1 1 8)  into  Eq.  ( 104),  we  obtain 


[  -  «7Ar(l  sin  Vn  +  ^(  -  ris^sin  Irff,  A', 


■f  [( -f  p®<a^)£,jv i4yy]sin  3i7^  AT,  +  ^  (— cr)\  +  p°(o^)A„  sin  i;„  AT, 


(119) 


and  we  note  that  the  solution  in  Eq.  (117)  with  Eq.  (118) 
satisfies  Eq.  (105)  identically  because  the  linear  terms  satisfy 
the  linear  Homogeneous  boundary  conditions  (71)  termwise 
and  the  nonlinear  term  satisfies  the  nonlinear  transcendental 
equation  (113)  Utilizing  the  fact  that  a  is  in  the  vicinity  of 
a>,v,  substituting  from  Eq.  (Ill)  into  Eq.  (119)  and  neglecting 
terms  of  order we  have 

[  -  «7^(1  +p]-¥pW]As  sin  Tf^  AT, 

=  (120) 
2h 


from  which,  with  the  aid  of  the  orthogonality  of  the 
sin  nirXy/V,  we  obtain 


{or  —  S}}f)A,^  1  =  0} 


-1)' 


iN-  11/2 


(121) 


and  we  do  not  bother  to  obtain  any  other  A„  because,  as 
already  noted,  we  are  interested  in  the  solution  only  when  oj 
is  in  the  vicinity  of  tujv  and  the  N  th  eigenmode  is  dominant  so 
that  all  other  eigenmodes  are  negligible.  Substituting  from 
Eqs.  (116)  and  (112b)  into  Eq.  (121),  we  obtain 

Ay  [o}^  —  a2^(l  +  ctAyA  *)] 

=  4{-ir-''^d,,{l/h)V„o}\  (122) 

where  Wy  is  given  in  Eqs.  (73)  and 


a  = 


—  9  y  ,  _ 


1)' 


fN-  1)/I 


"“I 


16  c  4/^ 

Equation  (122)  gives  the  nonlinear  relation  between  the  am¬ 
plitude  of  the  extensional  displacement  of  the  rod  and  the 
voltage  across  the  crystal.  Thus  for  an  <u  in  the  vicinity  of  <2;^ , 
the  solution  can  very  accurately  be  written 


Aysm  —  X, - _Ar,je*“', 


(124) 


where  must  satisfy  Eq.  (122).  Equations  (122)  and  (124) 
constitute  the  solution  of  the  steady-state  nonlinear  forced 
extensional  vibration  problem  for  the  rod  for  driving  fre¬ 
quencies  in  the  vicinity  of  o}y . 

Since  we  wish  to  obtain  the  nonlinear  relation  between 
the  amplitude  of  the  extensional  motion  of  the  rod  and 
the  driving  voltage  Vg  shown  in  the  circuit  in  Fig.  5,  we 
substitute  frmn  Eqs.  ( 1 24)  into  Eq.  (66),  which  is  then  substi¬ 
tuted  into  Eq.  (63)  to  obtain 


r  — 

4  =/<y4u>-^^F„ 


2A 

/ 


4. 

fJs-y.i 


(125) 


where  we  have  made  the  identification  4  =  —  /  since  the 
rod  is  passive'^  and 


Ay=Ay{-l] 


(AT  -  l)/2_ 


126) 


Application  of  KirchhofTs  voltage  equation  to  tl  rcuit 
shown  in  Fig.  5  yields 


+  4('l^g  +.1^2.)+  4,  =0,  ’27) 

where  Vg  is  the  generator  voltage,  Rg  and  R[_  are  tl 
ator  and  load  resistance,  respectively,  and  is  tht  '  ge 
across  the  crystal.  Substituting  from  Eqs.  (125)  and  (127)  into 
Eq.  (122),  we  obtain 


^^[o}^  ~al[\  -|-aA^AS)]^l  +  i<o  ^’^^^—{Rg  -|-/?^)j 

-  io)(Rg  +Ri,)4wj- efi {k i,  fAy  =  -  ^  F, , 

(128) 

and,  as  usual,  we  have  replaced  Oy  by  ^y  where 

0}y  =  0)y  -1-  io)y/2Qy,  (129) 


in  which  Qy  is  the  unloaded  quality  factor  of  the  resonator  in 
the  TV  th  mode.  Equation  ( 128)  is  the  nonlinear  equation  relat¬ 
ing  the  amplitude  of  the  extensional  displacement  to  the  gen¬ 
erator  voltage.  Since  Eq.  ( 128)  conteins  A  %  as  well  as  A  it  is 
not  simply  a  complex  cubic  equation.  In  order  to  facilitate 
the  solution  of  Eq.  ( 1 28),  we  multiply  Eq.  ( 1 28)  by  its  complex 
conjugate  and  after  some  algebraic  manipulation  obtain 

a\d^  -1-  -  2{d  +  Ca)a%  a®  -f  (1  -|-  a^)co*y  a4^] 

=  e\-g,  (130) 


where 

^=AyAX,  ,.g  =  VgV*,  (131) 

a  =  o}{Rg  +Rl  ){2wl€lj  /h),  +  o}%  a/Qy , 

C=  -  \<o\/Qy)  +  a{0^  -0}y)  —  b,  e  =  d^J/hv, 


V  =  /V  V/4^u^  b  =  Q>{Rg  +  R^)(Awleli/h  )(k  )Vv.  (132) 


Equation  ( 1 30)  is  a  real  cubic  algebraic  ecpiatipn  in  the  square 
of  the  extensional  vibratjon amplitude  Ay  A  which  may 
readily  be  sol ved  for  the  A  A  .  For  a  solution  A  y  A  J  to  be 
meaningful,  it  must  be  real  and  positive.  Consequently,  there 
is  at^  least  one  and  possibly  three  physically  meaningful 
Ay  A  X-  When  a  physically  meaningful  Ay  AX  has  been  de¬ 
termined  from  Eq.  (130),  it  may  be  substituted  in  Eq.  (128), 
which  then  becomes  a  linear  equation  in  Ay  that  may  readily 
be  solved  for  the  phase  if  desired.  When  an  Ay  is  known,  the 
current  4  through  the  circuit  may  be  evaluated  from  Eq. 
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(125).  This  result  can  be  used  as  a  check  for  the  y  determined 
from  intermodulation  measurements  for  any  desired  orien¬ 
tation  and  in  particular,  for  the  practical  5°  X  cut  of  quartz. 
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APPENDIX 

In  this  Appendix  we  derive  the  expression  in  Eq.  (49) 
with  Eq.  (50)  for  the  local  orthogonal  transformation  con¬ 
taining  terms  up  to  quadratic  in  the  mechanical  displace¬ 
ment  gradients.  To  this  end  we  first  note  that  Green’s  defor¬ 
mation  tensor  is  defined  by’® 

(Al) 

and  that  the  polar  decomposition  theorem  takes  the  form^' 

y,L=R>KVKL<  (A2) 

where 

^Kt4  ~  (A3) 

and,  of  course,  if  U  is  the  reciprocal  of  we  have 

RiK  =yu  uck-  (A4) 

We  now  substitute  from  Eqs.  (4)  and  (37a)  into  Eq.  (A1 ) 
to  obtain 

^KL  —  ^KL  +  +  ^i.K  ^i,L  •  (A5) 

Since  the  magnitude  of  small,  on  account  of  Eqs.  (A3) 

and  (A5)  U^m  must  differ  from  the  identity  8,.^  by  a  small 
(but  nonlinear)  quantity  thus 

Vkm  —  +  Ykm-  (A6) 

Substituting  from  Eq.  (A6)  into  Eq.  (A3)  and  employing  Eq. 
(A5),  we  obtain 

YkL  —  ^KL  +  i  ^i,K  ^I.L  ~  J  YkM  YmL  ■  (A7) 

Also,  since  small,  U  must  diffei  from  S^l  hy  a 

small  (but  nonlinear)  quantity  a^L  >  thus 

^  ML  —  ^ML  +  f^ML  •  (A8) 

Since 

VKMUiii.=8^^,  (A9) 

the  substitution  of  Eqs.  (A6)  and  (A8)  in  Eq.  (A9)  yields 

Okl  =  -Ykl-  Ykm  a  ml  ■  (A  10) 

Furthermore,  from  Eq.  (A4),  with  Eqs.  (4)  and  ( A8),  we  have 

Rik—^imRmk’  (All) 

where 

RmK  =  ^MK  +  ^M.K  +  ^LK’  (A  12) 

and  we  have  employed  the  translation  operator  dis¬ 

cussed  in  Sec.  I  in  obtaining  Eqs.  (All)  with  (A  1 2).  Equation 
(A  1 2)  shows  that  if  we  can  find  to  any  order  in  the  small 

quantities  u^xf  can  obtain  Rh,k  and,  of  course.  Rig  to 
that  order  in  u^x  ■ 

To  lowest  order,  from  Eqs.  ( A7)  and  (A  10),  successively. 


we  obtain 

oYkl  (A  13a) 

(Akl  —  ~  ^KL’  (A13b) 

the  substitution  of  the  latter  of  which  in  Eq.  (A  12)  yields  the 
usual  linear  relation  given  in  Eq.  (52).  To  second  order,  using 
(A13)  in  the  nonlinear  terms  in  Eqs.  (A7)  and  (AlO),  we  ob¬ 
tain 

iYkL  —  ^KL  +  J  ^M.K  ^MX  ~  J  ^KM  ^ML<  (A14a) 

\*^KL  —  ~  ^KL  ~~  i  ^M.K  ^MX  +  ^  ^KM  ^ML  ■  (A  14b) 
Substituting  from  Eq.  (A  13b)  into  the  nonlinear  terms  and 
(A  14b)  into  the  linear  terms  in  (A  12)  and  employing  Eqs. 
(37a)  and  (51),  we  obtain 

Rmk  —^mk  ^km>  (A15) 

where 

^KM  ~  ^KM  +  i  ^KL  ^LM  +  J('^MZ.  ^LK  ~  ^KL  ^ LM )' 

(A16) 

which,  respectively,  correspond  to  Eqs.  (49)  and  (50).  It  is  a 
simple  matter  to  show  that  the  orthogonality  relations 

RlK  RlM  —  ^KM<  RlK  RmK  —  ^LM  (A17) 

are  satisfied  to  second  order  in  the  u^^f. 
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'In  this  work  we  are  interested  in  obtaining  the  extensional  equation  for  the 
rod  essentially  in  the  limit  as  the  cross-section  collapses  to  a  line. 
"’Reference  5.  Chap.  7,  Sec.  1 . 

"R.  Bechmann,  “Elastic  and  Piezoelectric  Constants  of  Alpha-Quartz," 
Phys.  Rev.  110,  1060(1958). 

’’R.  N.  Thurston,  H.  J.  MeSkimin,  and  P.  Andreatch,  Jr.,  “Third-Order 
Elastic  Coefficients  of  Quartz,”  J.  Appl.  Phys.  37,  267  ( 1965). 

’  ’Reference  5,  Chap.  15,  Sec.  1. 

'"Only  terms  with  resonant  denominators  for  driving  frequencies  in  the  vi¬ 
cinity  of  the  resonant  frequencies  are  targe  and  need  be  retained  in  the 
nonlinear  products. 

’’Reference  2,  Sec.  II.  These  are  the  only  frequency  combinations  occurring 
at  the  lowest  orders  of  nonlinearity  that  are  in  the  vicinity  of  eo^. 

"Only  one  term  is  retained  because  only  the  largest  term  is  considered  to  be 
of  any  potential  importance.  If  more  are  thought  to  be  required,  they  may 
be  included  without  difficulty. 

’’The  procedure  is  a  straightforward  iterative  procedure  except  when  the 
functional  form  resulting  from  the  nonlinear  terms  is  identical  with  that 
satisfying  the  basic  linear  differential  equation,  at  which  point  the  coeffi¬ 
cients  associated  with  the  functional  form  due  to  the  nonlinearity  are 
treated  as  unknown  and  identical  with  the  coefficients  of  the  linear  solu- 
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tion.  This  is  the  reason  Eq.  (104)  is  a  homogeneous  equation  whereas  Eq.  '*It  should  be  recalled  that  the  last  term  in  Eq.  1 1 16)  is  the  large  one  for  the  . 

(99)  is  an  inhomogeneous  equation.  This  procedure  is  appropriate  because  case  S  =  I  only.  For  other  values  of  N,  values  ofM  have  to  be  consid- 

the  term  on  the  rhsof(  104)  actually  is  a  homogeneous  term  and  appears  to  ered,  perhaps  more  than  one.  However,  if  for  a  particular  orientation  k 

be  inhomogeneous  (known)  only  because  of  the  nature  of  the  formal  proce-  turns  out  to  be  very  small,  the  last  term  in  Eq.  (116)  is  negligible  and  the 

dure  employed  which  tends  to  be  misleading  whenever  this  type  of  coinci-  solution  holds  for  any  odd /V. 

dence  of  functions  occurs.  Truesdell  and  R.  A.  Toupin,  “The  Classical  Field  Theories,"  in  Ency-  -  - 

"A  similar  effect  was  inadvertently  neglected  in  Ref  3  and  shou'd  have  clopedia  of  Physics,  edited  by  S.  Fliigge  (Springer- Verlag,  Berlin,  I960),  ^  41 

been  included.  However,  the  quadratic  coefficient  in  Ref  3  for  the  AT-cut  Vol.  III/l,  Eq.  (26.2). 

turns  out  to  be  so  small  that  the  effect  has  a  negligible  influence  on  the  ^'Reference  20,  Eq.  (37.8).  The  Appendix  does  no  more  than  carry  out  to 
numerical  results.  second  order  what  the  reference  states  can  be  done. 


